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Abstract 

A simulation of an atomic force microscope operating in the constant amplitude dynamic mode 
is described. The implementation mimics the electronics of a real setup which includes a digital 
Phase Locked Loop (PLL). The PLL is not only used as a very sensitive frequency detector, but 
also to generate the time-dependent phase-shifted signal which drives the cantilever. The optimum 
adjustments of individual functional blocks and their joint performance in typical experiments are 
determined in details. Prior to testing the complete setup, the performances of the numerical 
PLL and of the amplitude controller were ascertained to be satisfactory compared to those of 
the real components. Attention is also focussed on the issue of apparent dissipation, that is 
of spurious variations in the driving amplitude caused by the non-linear interaction occurring 
between the tip and the surface and by the finite response times of the various controllers. To 
do so, an estimate of the minimum dissipated energy which is detectable by the instrument upon 
operating conditions is given. This allows to discuss the relevance of apparent dissipation which 
can be conditionally generated with the simulator in comparison to values reported experimentally. 
The analysis emphasizes that apparent dissipation can contribute to the measured dissipation 
up to 15% of the intrinsic dissipated energy of the cantilever, but can be made negligible when 
properly adjusting the controllers, the PLL gains and the scan speed. It is inferred that the 
experimental values of dissipation reported cannot only originate in apparent dissipation, which 
favors the hypothesis of "physical" channels of dissipation. 

PACS numbers: 07.79.Lh, 07.50.Ek, 46.40.Ff 

Keywords: virtual machine, non-contact AFM, dissipation, damping, apparent dissipation. Phase Locked 
Loop 
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I. INTRODUCTION 



Since almost a decade, non-contact atomic force microscopy (nc-AFM) has proven capable 
of yielding images showing contrasts down to atomic scale on metals, semiconductors, as well 
as insulating ionic crystals, with or without metallic or adsorbate overlayeraii^i^. Like other 
scanning force methods, the technique relies on a micro-fabricated tip grown at the end of 
a cantilever. However, unlike the widely used contact or the tapping modes, the cantilever 
deflection is neither static nor driven at constant frequency, but is driven at a frequency 
/o = a^/27r equal to its fundamental bending resonance frequency, slightly shifted by the 
tip-sample interaction. A sufficiently large oscillation amplitude prevents snap into contact. 
A quality factor exceeding 10^, readily achieved in UHV, together with frequency detection 
by demodulation provide unprecedented force sensitivity^"^. A phase-locked loop (PLL) is 
typically used for that purpose. Since /o varies with the tip-surface distance, it deviates from 
/o, the fundamental bending eigenfrequency of the free cantilever. Upon approaching the 
surface, the tip is first attracted, in particular by Van der Waals forces, which decrease /q. 
The negative frequency shift. A/ = /o — /o, varies rapidly with the minimum tip-distance 
d, usually as c?^" with n > 1.5, and then as exp{—d/X) a few angstroms above the surface, 
owing to short-range chemical and/or steric forces^. When A/ is used for distance control, 
contrasts down to the atomic scale can be achieved. Another specific feature of the nc-AFM 
technique is that the oscillation amplitude A is kept constant while approaching or scanning 
the surface at constant A/. 

Controlling the phase of the excitation so as to maintain it on resonance and to make the 
frequency matching a preset /o-value, as well as the driving amplitude so as to keep the tip 
oscillation amplitude constant, respectively, requires dedicated electronic components. Am- 
plitude control is usually achieved using a proportional integral controller (PIC), hereafter 
referred to as APIC, whereas phase and frequency control can be performed in two ways. 
In both cases the AC defiection signal of the cantilever is filtered, then phase-shifted and 
multiplied by the APIC output and by a suitable gain. The most common method consists 
in using a band-pass filtered defiection signal^"^. This is referred to as the self-excitation 
mode. The second method, extensively analyzed hereafter, consists in using the PLL to 
generate the time-dependent phase of the excitation signal. The PLL output is driven by 
the AC defiection signal and phase-locked to it, provided that the PLL settings are properly 
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adjusted. Then, the PLL continuously tracks the oscillator frequency /o with high precision. 
Moreover, the phase lag introduced by the PLL itself can be compensated. For reasons of 
clarity, this mode will be referred to as the PLL-excitation mode. The choice of the PLL as 
the excitation source has initially been motivated to take benefit of the noise reduction due 
to PLLs^. A further advantage is that the noise reduction does not only optimize the detec- 
tion of the frequency shift, but also the excitation signal. In both modes, the phase shifter is 
adjusted so as the phase lag ^ between the excitation and the tip oscillation equals — 7r/2 rad 
throughout an experiment. If all adjustments and controls were perfect, the oscillator would 
then always remain on resonance. 

The nc-AFM technique therefore requires the simultaneous operation of three controllers : 
PLL, APIC and distance controller, which keeps constant a given A/ while scanning the 
surface. Since the tip-surface interaction makes the dynamic of the oscillator non-linear, the 
combined action of those three controllers becomes complex. They can conditionally inter- 
play^i^ and therefore influence the dynamics of the system. Consider for instance the time 
the PLL spends to track /o is long compared to the time constant of the APIC. Then, the 
cantilever is no longer maintained at /o, but at a frequency slightly higher or lower. Conse- 
quently, the oscillation amplitude drops^'' and the APIC increases the excitation to correct 
the amplitude reduction. Such an apparent loss of energy, which can as well be interpreted 
as a damping increase of the cantilever, does not result of a dissipative process occurring 
between the tip and the surface, but is the consequence of the bad tracking of /q. So-called 
apparent dissipation (or apparent damping) remains under discussions in the nc-AFM com- 
munity, which hinders the quantitative interpretation of the experimental proofs of dissi- 
pative phenomena on the atomic scale over a wide variety of samples^'^'^'^'^'^'^ii^'i^'^. 
Thus, addressing the problem of apparent dissipation turns out to be mandatory but requires 
to understand the complex interplay between controllers as well as to analyze the system 
time constants. Although several models of physical dissipation, connected or not to the 
conservative tip-surface interaction have been proposed2ii2^i2^i2ii2^i2^i2Ii22ii2£i22 and reviewed^, 
the question of apparent dissipation in the self-excitation scheme has been addressed by two 
groups^'^'^. M. Gauthier et al. 7|] emphasize the interplay between the controllers and the 
conservative tip-sample interaction which, although weak, can significantly affect the damp- 
ing. They put in evidence resonance effects which can conditionally occur in damping images 
upon scan speed and APIC gains. G. Couturier et al.- address a similar problem numerically 



4 



and analytically. They show that the self-excited oscillator can be conditionally stable within 
a narrow domain of Kp and Ki gains of the APIC, but that consequent damping variations 
can as well be generated upon conservative force steps which change the borders of the sta- 
bility domain. The results mentioned above are valid for the self-excitation mode, but the 
question of apparent dissipation remains open regarding the PLL-excitation mode. However 
recently, J. Polesel-Maris and S. Gauthier [ssj have proposed a virtual dynamic AFM based 
on the PLL-excitation scheme. Their work is targeted at images calculations including real- 
istic force fields obtained from molecular dynamics calculations^i. Their conclusions stress 
the contribution of the scanning speed and of the experimental noise to images distortion 
but do not address the potential contribution of the PLL upon operating conditions. 

The goal of the present work is two-fold : 1- providing a detailed description of the 
PLL-excitation based electronics of a home-built AFM used in our laboratory; 2- assessing 
the contribution of the various controllers to the dissipation signal and in particular the 
contribution of the PLL. The paper is organized as follows. In section [Tll an overview of 
the chart of the microscope and of the attached electronics (c/. figH]), is given in terms 
of blocks, namely ; oscillator and optical detection (block 1), RMS-to-DC converter (block 
2), amplitude controller (block 3), PLL (block 4), phase shifter (block 5) and tip-surface 
distance controller (block 6). In section UTTt the detailed description of the numerical scheme 
used to perform the calculations is given on the base of coupled integro-differential equations 
ruling each block. Section IIVI provides an estimate of the minimum detectable dissipation 
by the instrument with the goal to assess the relevance, compared to experimental results, 
of the apparent dissipation which can be conditionally generated numerically. Section |V] 
reports the results. In the first part, the simulation is validated by comparing a numerical 
A/ vs. distance curve to the analytic expression of the A/ due to Morse and Van der Waals 
interactions which does not take into account the finite response of the various controllers. 
Then the dynamic properties of the numerical PLL and APIC upon gains are compared to 
those of the real components. Section IV CI gives some examples on how apparent dissipa- 
tion can be produced upon working conditions of the PLL. Section IV Dl finally shows scan 
lines computed while varying PLL gains, scan speed and APIC gains. A discussion and a 
conclusion end the article. 
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II. OVERVIEW 



A. Description 

The electronics consists of analog and digital (12 bits) circuits which are described by 
six interconnected main blocks operating at various sampling frequencies (/s). The highest 
sampling frequency among the digital blocks is the PLL one, fg^ = 20 MHz. The PLL 
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electronics has initially been developped by Ch.Loppacher 

Block 1 represents the detected oscillating tip motion coupled to the sample surface. In 
the simulation, the block is described by an equivalent analog circuit. More generally, all the 
analog parts of the electronics are described in the simulation using a larger sampling fre- 
quency compared to /s^, namely = 400 MHz. This is motivated by the ultra-high vacuum 
environment within which the microscope is placed, thus resulting in a high quality factor 
of the cantilever, typically Q = 30000 at room temperature. Besides, nc-AFM cantilevers 
have typical fundamental eigenfrequencies /o — 150 kHz. The chosen sampling frequency 
should therefore insure a proper integration of the differential equations with an error weak 
enough. The signal of the oscillating cantilever motion goes into a band pass filter which 
cuts-off its low and high frequencies components. The bandwidth of the filter is typically 
60 kHz, centered on the resonance frequency of the cantilever. Despite the filter has been 
implemented in the simulation, no noise has been considered, so far. The signal is then sent 
to other blocks depicting the interconnected parts of two boards, namely an analog/digital 
one, the "PLL board", and a fully digital one which integrates a Digital Signal Processor 
(DSP), the "DSP board". The boards share data via a "communication bus" operating at 

= 10 kHz, the lowest frequency of the digital electronics. 

Block 2 stands for the lone analog part of the PLL board (/^ = fsz)- It consists of a RMS- 
to-DC converter. The block output is the rms value of the oscillations amplitude, y4rms(^)- 
^rms(^) is provided to block 3, one of the two PlCs implemented on the DSP (/« = fg.^)- 
When operating in the nc-AFM mode, the block output is the DC value of the driving 
amplitude which maintains constant the reference value of the oscillations amplitude, A^^^. 
This is why it is referred to as the amplitude controller, APIC. For technical reasons due to 
the chips, the signal is saturated between and 10 V. 

The dashed line in figiH depicts the border between analog and digital circuits in the 
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PLL board. The digital PLL, block 4 (/^ = /^J, consists of three sub-blocks : a Phase 
Detector (PD), a Numerical Controlled Oscillator (NCO) and a filtering stage consisting of 
a decimation filter and a Finite Impulse Response (FIR) low pass filter in series. The PLL 
receives the signal of the oscillation divided by v4i.jns(^) plus an external parameter : the 
"center frequency", /cent = ^^cent/27r. /cent specifies the frequency to which the input signal 
has to be compared to for the demodulation frequency stage. This point is particularly 
addressed in section IIIIDI The NCO generates the digital sin and cos waveforms of the 
time-dependent phase, fnco(t) + iPpn(t), ideally identical to the one of the input signal. fpn{t) 
is correlated to the error which is potentially produced while the frequency demodulation, 
upon operating conditions. The sin and cos waveforms are then sent to a digital phase 
shifter, block 5 (/^ = /sj which shifts the incoming phase (pncoit) + v^pii(t) by a constant 
amount, ipps, set by the user. Since the cantilever is usually driven at /o, v^ps is adjusted to 
make that condition fulfilled^, namely : 

^nco{t) + (Pp\\(t) + ipps = uJot, (1) 

Indeed, the PLL produces the phase locked to the input, that is ipncoit) +V^pii(''^) — uJ^t — n/2. 
If it optimally operates, v^pii(t) — 0. ipps has therefore to be set equal to +7r/2 to maintain 
the excitation at the resonance frequency prior to starting the experiments. Consequently, 
ip = — 7r/2 rad. The block output, sin [(Pncoit) + fpiiit) + V^ps], is converted into an analog 
signal and then multiplied by the APIC output, thus generating the full AC excitation 
applied to the piezoelectric actuator to drive the cantilever. 

Block 6 is the second PIC of the DSP (/^ = fs-^)- It controls the tip-surface distance to 
maintain constant either a given value of the frequency shift, or a given value of the driving 
amplitude while performing a scan line (switch 3 set to location "a" or "b", respectively 
in figil]). The output is the so-called "topography" signal. The block is referred to as the 
distance controller, DPIC. 

Finally, a digital lock-in amplifier detects the phase lag, (p, between the excitation signal 
provided to the oscillator and the oscillating cantilever motion. 
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B. Time considerations 



Analog and digital data are properly transformed by Analog-to-Digital and Digital-to- 
Analog Converters (ADC and DAC, respectively). In the electronics, ADCl is an AD9042 
{cf. figH]) with a nominal sampling rate of 41 x 10^ samples per second^''. This ensures the 
analog signal is sampled quick enough and properly operated by the PLL at /s^. This ADC 
is therefore not described in the simulation. ADC2 {ADS 7805) has a nominal frequency of 
100 kHz [37]. The signal is transmitted to the communication bus, the bandwidth of which 
is ten times smaller. Its role is therefore as well supposed to be negligible. The code is 
implemented assuming that the RMS-to-DC output signal is provided to the communication 
bus operating at /^g. DACl {AD 668) is a 12 bits ultrahigh speed converter. It receives the 
digital waveform coming from the PS. Indeed, it must be fast enough to provide a proper 
analog signal to hold the excitation. Its nominal reference bandwidth is 15 MHz 37|. To 
make the code implementation easier, the DAC has not been implemented neither. Thus, 
it is assumed that the PS signal directly provides the signal at /^^ to perform the analog 
multiplication, itself processed at fg,^ due to the APIC output. The others DACs have all 
nominal bandwidths much larger than the communication bus one and are also assumed to 
play negligible roles. 



III. NUMERICAL SCHEME 



A. Block 1: oscillator and optical detection 

The block mimics the photodiodes acquiring the signal of the motion of the oscillating 
cantilever. The equation describing its behavior is given by the differential equation of the 
harmonic oscillator : 

m + + u;lz{t) = ulEUt) + ^^^^ (2) 

ujq = 271 fo, Q, kc stand for the angular resonance frequency, quality factor and cantilever 
stiffness of the free oscillator, respectively. z{t), Hexc(i) and -Fint(t) are the instantaneous 
location of the tip, excitation signal driving the cantilever and the interaction force acting 
between the tip and the surface, respectively. The equation is solved with a modified Verlet 
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algorithm, so-called leapfrog algorithnt^S,^ using a time step At^^ = l//s2 = 5 ns. In the 

foUowings, the time will be denoted by its discrete notation : t ^ ti = i x Ats2- 

The instantaneous value of the driving amplitude Sexc(^i) (units : m) can be written as : 

(units : m.V~^) represents the linear transfer function of the piezoelectric actuator driving 
the cantilever. A^^citi) (units : V) is the APIC output (c/. section [Til Cp . It is proportional 
to the damping signal according to : 

KsAUU) = (4) 

T(ti) and Aq (units : s~^ and m, respectively) being the damping signal and oscillations 
amplitude of the cantilever when driven at /o, respectively. When the cantilever is externally 
driven and if no interaction occurs, Aexc(^i) can be written as a function of Aq and of the 
quality factor of the cantilever : 

i^3^xc,0 = ^ (5) 

Then the damping of the free cantilever equals : 

To ^ ^ (6) 

In nc-AFM, the dissipation is commonly expressed in terms of dissipated energy per oscil- 
lation cycle, Edg. For a cantilever with a high quality factor oscillating with an amplitude 
Ao : 

7ik,Al nk.A^To 
Edo[Ao) = = ■ [7) 

In UHV and at room temperature, Q = 30000. Besides, nc-AFM commercial cantilevers 
have typical stiffnesses^^ fc^ ~ 40 N.m~^. Considering Aq = 10 nm, the intrinsic dissipated 
energy per cycle of the cantilever is then Ed^ — 2.6 eV/cycle. 

In equation [3l Zps{ti) is the AC part of the excitation signal [of. section IIIIEI) . It is 
provided by the PS when the PLL is engaged. When the steady state is reached, e.g. 
ti ^ tstcady — '^Q / fo^ the block output is : 
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Kiz{ti) = KiA{U) sin [ujU + viU)] (8) 

Ki (V.m~^) depicts the transfer function of the photodiodes which is assumed to be hnear 
within the bandwidth (3 MHz in the real setup). If the damping is kept constant, the amph- 
tude and the phase, A(ti) and ip(ti) respectively, are supposed to be constant as well. This 
is no longer true once the various controllers are engaged, therefore their time dependence 
is explicitly preserved. 

In equation [21 the interaction force -Fint('") = ~<9r^nt('") is derived from a conservative 
potential consisting of two components : a long-range part, depicted by a Van der Waals 
term defined between a sphere and a half-plane and a short-range part, prevailing at closer 
distances, depicted by a Morse potential : 



or 



_ r — rc _ 2(r — rc) 

2e ^ — e ^ 



(9) 



H and R are the Hamaker constant of the tip-vacuum-surface interface and tip's radius, 
respectively. Uq, Tc and A are the depth, equilibrium position and range of the Morse 
potential. The instantaneous tip-surface separation is r(tj) = -D(tj) — z{ti), where D{ti) is 
the distance between the surface location and the cantilever position at rest. So far, neither 
elastic deformation of the sample and tip, nor dissipative interaction have been considered. 

The signal Kiz{ti) then gets into the band pass filter (EPF), the central frequency of 
which, /c = Wc/Stt, equals the resonance frequency of the cantilever, /o, with a bandwidth 
Bw — 60 kHz. The output, Zbpf(^j) (units : V), is ruled by : 

i;bpf(t) + 27rSw/ibpf(t) + c^c^bpf(t) = 2'KBwKiz{t) (10) 
Zhpi{ti) is then provided to the RMS-to-DC converter of the PLL board. 



B. Block 2: RMS-to-DC converter 



The converter is the only analog part of the PLL board. The related differential equation 
is integrated at /^j. The chip (AD734) computes the square root of the squared value of the 
incoming signal, preliminary filtered by a first-order low pass filter, the cut-off frequency of 



10 



which is fco = 400 Hz. The output is the amphtude (DC value) of the oscillation, Aj-^siti) 
(units : V) : 



(11) 



Vs(ti) being the output of the first-order low pass filter : 



TrmsVs{t)+Vs{t)=zljti), 



(12) 



with Trras = l/(27r/co) ~ 400 /iS. 

Zhpf{ti) is then divided by A,.^s{ti) in order to normalize the amplitude of the waveform. 
The signal thus normalized is sent to the ADCl to be operated by the digital PLL. 

C. Block 3: amplitude controller 

The block represents a digital PI controller implemented in the DSP board. The con- 
troller receives the RMS-to-DC output signal via the communication bus. Since the bus 
operates at fs^ = 10 kHz, the time step used to solve the related differential equation is 
Ats3 = l//s3 = 100 fis. Besides A^-^s, the controller receives three external parameters : 
the proportional and integral gains, K^'^ and K^'^ respectively (units : dimensionless and 
s~^, respectively), and the reference amplitude expected to be kept constant as soon as the 
controller is engaged (switch 1 set to location "b" in figH]), A'^^ (units : V). The block output 
is the DC value of the excitation, previously referred to as AcxdU) {cf. equj3]) : 



Engaging the APIC makes the nc-AFM mode effective. This requires the PLL-excitation 
mode (block 4, cf. section [Tll Dl) to be already engaged. If operating at /o, then Aj-^s/Ki 
equals the resonance amplitude, Aq. A^xc is then minimal. 



^exc(^i) —Kp^ [-^o'^ ~ A^jnsiti)~\ 




(13) 
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D. Block 4: PLL 



Before starting this section, note that some of the elements detailed hereafter are adapted 
from the book by R.Best 9|. The digital PLL consists of three sub-units : a Phase Detector, 
a decimation filter and a FIR low pass filter in series and a NCO. The block operates at 
fsi = 20 MHz, with the related time step Atg^ = 1//^^ = 50 ns. In the electronics, various 
FIR low pass filters have been implemented upon the desired sensitivity in the frequency 
detection, among which a 19**^ order filter with a 3 kHz cut-off frequency and a 45**^ order 
filter with a 500 Hz cut-off frequency. Both of them can be used in the simulation. 



1. Phase detector 



The PD is analogous to a multiplier regarding the two input signals : the BPF output 
divided by Arms and the cos waveform coming out of the NCO (c/. figil]). Their product is 
multiplied by a further gain, Kd (units : V) converting the dimensionless signal into volts 
to be operated by the FIR low pass filter. The instantaneous block output is referred to as 
KdZeiU) : 

KdzM = Kd^^^ cos (^nco(tO) (14) 

2. Filtering stage 

Assume that uJ^iti) and tUncol^i) are the instantaneous angular frequencies of the can- 
tilever and of the signal generated by the NCO, respectively. KfiZ(.{t.i) consists of a high 
frequency component : uJ^iti) + oj^coiti) and a low frequency one : 5uj{ti) = u5o(ti) — 
i^nco(^i)- The FIR low pass filter cuts off the high frequency component and produces 
Uf(ti) oc sm{6u!(ti)ti} oc [5c<j(i(:j)] x tj, which can be referred to as an error signal of the 
PLL. Indeed, when the PLL optimally operates, tUncol^i) almost perfectly matches uJo{ti). 
The instantaneous value Uf{ti) can therefore be interpreted as a correction term in the PLL 
cycle. 

Before being operated by the FIR low pass filter, the signal is processed by the decimation 
filter. The filter averages KdZe{ti) over Nds PLL cycles upon the FIR low pass filter cut-off 
frequency. For instance, A^^^ = 400 for the 3 kHz low pass filter. The updating rate of the 
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FIR low pass filter is therefore fs^/Nds = 50 kHz. The digital data are averaged over those 
Nds cycles. The average value is fed at the first entry of a buffer B consisting of A^fir entries. 
The entries of the buffer are then all shifted by one into the buffer. At a given moment in 
time, tj, Uf{ti) is given by the following algorithm : 



A^fir is the order of the FIR low pass filter (A^gi. <^ A^^) and is the k*^ coefficient of the 
FIR low pass filter. Once the buffer is transmitted, it is initialized and filled again. Finally, 
Uf{ti) is multiplied by a further gain, Kq, which depicts the linear conversion of the signal 
from volts to rad.s"^ (units: rad.V~^.s~^) and provided to the NCO. 

3. Numerical Controlled Oscillator 

We first assume that the frequency tracker of the PLL is disengaged (switch 2 set to 
location "a" in fig HI). Its role is carefully addressed in section UlI D 51 The NCO adds the 
instantaneous angular frequency KQUf{ti) to an external input, the center angular frequency 
of the PLL, 

ujcent- lucent IS fixed equal to the angular resonance frequency of the free cantilever 
ujq, prior to starting the experiments. The signal is then integrated, which produces the 
related phase, v^nco(^i), locked to the one of the cantilever : 



fc=pii 

tpu being the moment when the PLL is engaged. Obviously, the PLL has to be engaged once 
the oscillator has reached its steady state and before the APIC. 

4- Frequency demodulation 

When the tip is located far from the surface, ujq = uq. Once approached close enough 
from it, uJq starts decreasing. Meanwhile, the NCO produces uj^coiti) = uJcent + KoUfiti), 
as mentioned above. When the frequency tracker is disengaged, cUcent is kept constant and 
matches the resonance frequency of the free cantilever, Ucent = "-^o- Therefore KQUf{ti) is 




(15) 




(16) 
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nothing but the instantaneous frequency shift (actually 27r x A/) of the tip interacting with 
the surface. In other words : 



Wcent + KoUf{ti) -UJo = 27lAf{ti) (17) 

Kq is a key parameter of the PLL. It sets its capability to get locked to the input signal and 
in turn it sets its stability. R.Best defines Kq from the locking range Aui of the PLL, e.g. the 
frequency gap with respect to the center frequency the PLL can detect remaining locked^. 
On the hardware level, the control signal Uf{t) is limited to a range which is smaller than 
the supply voltages, usually ±5V. Assuming Uf^^^ and Uf^ be the minimum and maximum 
values allowed for Uf, Best defines Kq as : 

Ko = (18) 

Therefore Kq is related to the maximum frequency shift detectable per volt within the 
detection range of the low pass filter. Practically, the value of Kq is not accessible a priori. 
It's easier to set the locking range Aui. For an oscillation at /o = 150 kHz, frequency shifts 
of about a few hundreds of hertz are typically expected'^'^. We can therefore choose the 
3 kHz FIR low pass filter to insure a proper detection of A/, which sets the locking range to 
AuJi = 27C X 6000 rad.s^^. The maximal value of Kq expected is then ^ 11000 rad.V^^.s^^, 
which is an excellent estimate as detailed in section |Vl 



5. Frequency tracker 

The frequency tracker is a specific feature of our digital PLL. When engaged (switch 2 
set to location "b" in figil]), the center frequency is continuously updated by the FIR low 
pass filter output : 



a;cent(^i) = ^cent(ii-l) + KQUf{ti) (19) 

The updating frequency is 2.5 kHz. The frequency tracker has been implemented in order 
to compensate the fact that the frequency demodulation was performed via the lone propor- 
tional gain Kq. Thus, as mentioned before, KQUf{t.j) can be interpreted as the error signal 
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produced in the frequency detection compared to cUcent- Consequently, this error is also inte- 
grated by the NCO, which leads to an additional phase lag added to ip^co at each PLL cycle 
and previously referred to as v^pu- v^pu per PLL cycle can approximately be estimated to : 

^Vpw = V5pii(i^j+i) - V^piK^i) ~ KoUfiU) X Ats^ (20) 

ippw would be zero if no frequency shift occurred, which is the case in most of the applications 
using PLLs. But while approaching, Af decreases continuously, therefore so does ippw. On 
the contrary, when the frequency tracker is engaged, cUcent is continuously updated. The 
error in the frequency detection drops to zero. More exactly, it is equal to the difference 
between two consecutive values of c^cent '■ e — t^centl^i) — ^ccnt(^i-i)) but is necessarily small 
and so is (pp\\. 

To assess how sensitive to frequency changes the phase is, the following experiment is 
carried out. A 150 kHz sinusoidal waveform is generated by means of a function generator 
and sent to the real PLL. The frequency is then slowly detuned from —150 Hz up to +150 Hz. 
The phase lag between input and output waveforms, ippw, is recorded with a lock-in amplifier 
(Perkin Elmer 7280) and reported as a function of the detuning. The PLL center frequency 
is fixed to /cent = 150 kHz. The experiment is repeated the frequency tracker being engaged 
and disengaged. Two amplitudes of the PLL input waveform are used. In this experiment, 
the input waveform stands for the oscillatory motion of the cantilever and the tuning for the 
shift occurring when the tip is approached towards the surface upon attractive or repulsive 
forces. The results are reported in figEl When the tracker is disengaged, the maximum 
detuning corresponds to a phase lag of ±80 degrees, which means that the cantilever would 
then be driven off resonance severely. On the opposite, when engaged, the phase lag reduces 
(inset) to ±0.05 degree. 

This feature has no consequence when the PLL is only used as a frequency demodulator 
like in the self-excitation mode. On the opposite in the PLL-excitation scheme, this point is 
crucial since the PLL is produces the excitation signal. Therefore particular attention has 
to be paid on the way it is produced. If it is abnormally phase shifted, then the oscillation 
amplitude drops and consequently apparent dissipation is generated, as shown in section |Vl 
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E. Block 5: phase shifter 



The PS receives the sin and cos waveforms generated by the NCO. A further input to the 
block is the phase lag, (pps, fixed prior to starting the experiments to make the cantilever 
oscillating at /q. The PS digitally computes : 

Zps{ti) = sin [iPncoiU) + ippw] X COS (V3ps) 

+ COS [(pncoiU) + Vpw] X sin (V?ps) 

= sin [(fncoiU) + ippU + V5ps] (21) 

When the system is being operated in the PLL-excitation mode, Zps{t) is converted into an 
analog signal by the DACl and multiplied by the APIC output. 

F. Block 6: distance controller 

The distance controller is the second digital PI controller implemented in the DSP oper- 
ating at /sg. The block gets the setpoint value of the signal (A/ or damping) onto which 
the control of the tip-sample distance is performed and the proportional and integral gains, 
Kp'^ and K'^'^, respectively. Here, let's assume that the reference signal is the frequency 
shift, as depicted in figUl We have arbitrarily chosen not to describe the transfer function 
of the z-piezo drive. Therefore K^"^ and Kf'^ have natural units (nm.Hz~^ and nm.Hz~^.s~^, 
respectively). The controller is described by : 

D{U) =D(tdc) + K^' [A/set - AfiU)] 

i 

+ Y,Kf^[AU-Af{tk)]At,,, (22) 

A;>dc 

D(tdc) being the tip-surface distance when the DPIC is engaged. 

G. Lock-in amplifier 

The description of the lock-in amplifier implemented in the simulation does not depict 
the detailed operational mode of the real lock-in which is used to monitor the phase shift 
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of the oscillator (Perkin Elmer 7280). Its purpose is to provide an easy way to estimate 
the phase shift between the excitation and the oscillation. The calculation of the phase is 
performed at 2.5 kHz. The buffer used to extract the phase therefore consists of niock-in = 
/s^/2.5 kHz = 8000 samples. The numerical code used to describe it is : 

taii{if{ti)) = 

El=i-n,„,k-in ^bpf(tfc) X COs[(^nco(4) + ^^pU + <^ps] 

H. Code implementation 

The numerical code has been implemented with LabView^^ 6.1, supplied by National 
Instruments"""^. It consists of a user interface where all the parameters are tunable at run- 
time, like during a real experiment. The couple of integro-differential equations [21 [3l [TOl 
[TTl W2i [131 [m [13 [13 [211 and [22] are integrated at their respective sampling frequencies. 
The monitored signals are the oscillation amplitude Arms (equ JTTl) . the frequency shift A/ 
(equJTTI). the phase ip (equl23!) and the relative damping T/TQ — l=QK3Aexc/Ao — ^, deduced 
from the APIC output (equH]). The connection to the dissipated energy per cycle is given 
by equation [71 that is P/Pq — 1 = Ed/ Ed^ — 1. 

IV. APPARENT DISSIPATION VS. MINIMUM DISSIPATION 

Addressing the question of apparent dissipation requires to estimate the minimum dis- 
sipation which is detectable by the instrument upon operating conditions. Beyond the 
specificities of the PLL- or self-excitation modes, important parameters like quality factor 
Q, temperature and bandwidth of the measurement must be considered. 

We here focus on the minimum dissipated energy, 5Ed, due to thermal fluctuations of 
the cantilever when it oscillates close to a surface. Thermal driving forces are connected 
to the energy dissipation by the Q factor of the cantilever. The thermal kicks introduce 
fluctuations of amplitude and phase and therefore fluctuations of the energy dissipation. 
This is true for a free cantilever, but the contribution of the thermal noise is expected to 
be even more pronounced when the tip is close to the surface. Then, the fluctuations of the 
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interaction force SFi^t have a strong influence on the nonhnear dynamics of the cantilever, 
in particular when the tip is at distances involving short-range forces where the nonlinearity 
is more pronounced. 



n 

The instrumental noise (c/. Ch.2 in refs.[l| and [33|), essentially due to electronic com- 
ponents, is not considered and we further assume that the electronic blocks (RMS-to-DC, 
PI controllers, PLL) operate perfectly. Doing so, SE^ is under-estimated but the framework 
of this section is to provide a ground value to be compared to the values obtained with the 
simulation. 

A. Connection between 6Efi and SF-mt 

The fluctuation of the dissipated energy per cycle can be connected to the fluctuation of 
damping 6T, via equJT] : 

5 Ed = TTkcAl— (24) 

UJq 

Besides, because Aexc = Aq/Q = AoTq/uq = Fexc/kc on resonance and because the tip- 
sample interaction force F[^t can be treated, to flrst order^-, on the same level as -Fexc, a 
fluctuation of -Pint should produce, a fluctuation of damping : 



_ (5Acxc _ SFir,t 

Consequently : 



(25) 



5Ed = TiAoSF^^t (26) 



B. Estimate of 8F\ 



int 



For large oscillation amplitudes (that is larger than the minimum tip-surface dis- 
tance, a few angstroms), Fint is connected to the so-called normalized frequency shift^, 
7 = AfkcA^^'^/fo, via the equation^ (c/. also Ch.l6 in ref. pj]) : 



7{r) ~ 0.43v/V^nt(r)Fint(r), (27) 
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where Vint{r) and Fi^t{i^) are the interaction potential and force, respectively, between the 
tip and the sample at a location r. The fluctuation in the relative frequency shift 5A///o = 
Sf / fo, that is the cantilever frequency noise, due to a fluctuation of Fint is then given by : 



^/.^ Mrl,^,. (28) 



C. Estimate of 6f/fo 



Y. Martin et al. 



in ref. 



44 , T.R. Albrecht et al. 



. fl,H. 



Diirig et al. |45| and F.J. Giessibl (Ch.2 



ij) have calculated the thermal limit of the frequency noise in frequency-modulation 



technique over a measurement bandwidth B. It is given by 



5f / 2kBTB 



(29) 



Therefore, the dissipated energy due to thermal fluctuations of the cantilever close to the 
surface can be estimated to : 



,^,.4.6.?MWiM (30) 

V 7r/oQMnt(r) 

The measurement bandwidth B can be estimated out of the following considerations. As 
mentioned by F.J. Giessibl (c/. Ch.2 in ref.|l|), i? is a function of the scan speed Vg and the 
distance Qq between the features which need to be resolved : 

B = - (31) 

For UHV investigations, Oq is of about one atomic lattice constant, that is a few angstroms. 
At room temperature, due to thermal drift, atomic scale images are usually recorded at scan 
speeds of about 6 lines (3 forwards plus 3 backwards) per second. Let's consider for instance 
a moderate resolution of 6 pixels per atomic period. Then, a line consisting of 256 pixels 
should be acquired with a bandwidth B = 6 x 256/6 = 256 Hz. 

Table [H gives some estimates of the relative dissipated energy due to thermal fluctuations 
of the cantilever dEd/E^Q close to the surface in the short-range or pure Van der Waals 
regimes at various temperatures and for various quality factors. In UHV at room tem- 
perature, our experimental conditions, the minimum dissipated energy which is detectable 
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Q Interaction E^^ SE^ dEd/E^g 

regime (cV/cydc) (oV/cyclc) 



5000 (298°K) 


VdW + short-range 


7.69 


0.177 




2.3% 




VdW only 


7.69 


0.141 




1.8% 


30000 (298°K) 


VdW -|- short-range 


1.28 


7.25 X 10" 


-2 


5.7% 




VdW only 


1.28 


5.78 X 10- 


-2 


4.5% 


500000 (4°K) 


VdW -f short-range 


0.077 


2.06 X 10" 


-3 


2.7% 




VdW only 


0.077 


1.64 X 10" 


-3 


2.1% 



TABLE I: Dissipated energy of the free cantilever E^q (equj?!) and dissipated energy due to thermal 
fluctuations of the cantilever close to the surface SE^ (equ l30]) for various quality factors and 
temperatures when Van der Waals plus short-range (equED or pure Van der Waals forces (similar 
equation, with Uq = 0) are considered. The cantilever parameters are Aq = 7 nm, /o = 150 kHz, 
kc = 40 N.m~^ and B = 260 Hz. The parameters of the interaction potential have been taken from 



ref. 3 : 



H = 1.865 X 10"^^ J, i? = 5 nm, Uq = 3.641 x 10"^^ J, A = 1.2 A, and rc = 2.357 A. 6Ed 
has been estimated at a distance r for which the two interaction regimes are clearly distinct (c/. 
figUa)), r = 5 1. 

corresponds to 5% of the intrinsic dissipated energy of the free cantilever. This corresponds 
to about 150 meV/cycle with typical conditions for UHV investigations carried out at room 
temperature (c/. equH and discussion below). Besides, as mentioned before, this value is 
underestimated. A straightforward consequence is that the strength of apparent dissipation 
should overcome this limit to be relevant. With a moderate quality factor in the Van der 
Waals regime like in high vacuum for instance, the limit drops by almost a factor 3 (1.8%). 
Thus, apparent dissipation effects might occur more easily under these conditions'^. At low 
temperatures, in the short-range regime the ratio is 2.67%. However, this value is likely 
still too high because then, the thermal drift being drastically reduced, the measurement 
bandwidth can be lowered and apparent dissipation more likely to be measured. 
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V. RESULTS 



A. Validation of the numerical setup 

Frequency shift vs. distance curves obtained from the simulation have first been com- 
pared to the analytic expression of A/ due to Van der Waals and Morse potentials (c/. 
appendix, section A). The results are shown in fig|3]^a). The parameters chosen to perform 
the simulation are consistent with typical parameters used during experiments performed 
in UHV. The parameters of the interaction potential have been taken from ref.j47|]. They 
are representative of the interaction between a silicon tip and a silicon(lll) facet. An ex- 
cellent agreement is observed between numerical and analytic curves along the attractive 
and repulsive parts of the interaction potential, thus validating the numerical scheme. The 
parameters used to perform the calculation are given in the caption. Let's also notice that 
the frequency tracker was engaged. In figs. Mjo), (c) and (d), the variations of ip, A^ms 
and relative damping, respectively are reported vs. the tip-surface separation. Phase and 
amplitude remain almost constant while approaching, within, however, deviations limited 
to 0.3% compared to —90 degrees and A^^^ = 7 nm, respectively. In the repulsive part of 
the potential, steep phase changes occur, but the amplitude does not dramatically drops, at 
least up to A/ = +100 Hz. Consequently, the relative damping remains constant. 

B. Numerical vs. real setups 

1. PLL dynamics 

The dynamic behaviors of real and simulated PLLs have then been compared. The 
experiment consists in locking the PLL onto a 150 kHz sinusoidal waveform according to 
the same procedure than in section HHP 51 The 3 kHz FIR low pass filter is used. At a 
certain moment, a frequency step of +10 Hz is applied to the center frequency, resulting in 
a shift of —10 Hz {uq + 27rA/ = cjcont)- The step response is recorded for various values of 
the so-called loop gain (real PLL) and various values of K^Kd (simulation). The variations 
of A/ vs. time are fitted with simple decaying exponential functions, the characteristic time 
of which stands for the locking time of the PLL. The results are reported in figs|l](a) and 
(b). A rather long locking time is noticed for low values of the gains whereas the PLLs lock 
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faster when the gains become larger. For the latter case, the PLLs can operate up to the 
limit of the locking range as shown by the oscillations. 

The locking time deduced from each fit is plotted as a function of the gains of both PLLs. 
In order to make the curves comparable, the loop gain must be rescaled by an arbitrary 
constant which depends on the electronics. The best agreement between the curves was 
achieved with 91 x 10^ (c/. fig|5]). A single master curve of the PLL dynamics can thus 
be extracted. The rather good agreement between the two curves provides evidence that 
the simulation reasonably describes the real component, at least within the locking range. 
For values of the gains up to K^Kd = 6000 rad.s~^, the PLL is stable and able to track 
frequency changes within the locking range around 150 kHz. Above 6000 rad.s^^, the PLL 
introduces overshoot in the output waveform while attempting to lock the input signal. 
For higher gains, the PLL is not able to properly track the input signal, even though its 
frequency is within the locking range. The border is reached for Kq > 10^ rad.V~^.s~^, in 
good agreement with the value expected from R. Best's criterion (c/. discussion in section 
IIIID4I) . The arrow in figJS] indicates the usual loop gain value which is chosen to perform 
the experiments using the 3 kHz low pass filter, corresponding to K^Kd = 5000 rad.s"^. The 
related locking time of the PLL is then ~ 0.35 ms. For those values of gains, the locking 
range is about ±400 Hz. 

2. APIC dynamics 

In order to extract a typical time constant of the component, similar experiments have 
been carried out with the APICs. The cantilever remaining far from the surface, a step is 
applied to the setpoint amplitude A^^ resulting in an abrupt change in Aj-ms upon gains. 
The results are reported in figl6](a, real setup) and (b, simulated setup). The curves exhibit 
over- (no overshoot at all) under- (oscillating behavior) or critically damped (single over- 
shoot) behaviors upon chosen gains. So as to extract the APIC response time, we focus 
at curves which exhibit a single time constant, that is curves for which a weak overcriti- 
cally or a critically damped response is observed (c/. insets in figsj6t^a) and (b)). This is 
motivated by the controller response which is then the fastest, while preserving an overall 
stable behavior. The changes in Arms are fitted with decaying exponentials functions and 
the related characteristic time is extracted. The variation of the so-called response time of 
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the controller (tresp) vs. gains is reported in fig|71 The restriction to curves exhibiting a 
single time constant is similar to restricting the analysis to a single gain of the controller—. 
Thus, a single master curve which describes the dynamics of both APICs can be extracted 
as well. In figl?! the Kp gain of the real controller has been rescaled to make it matching 
Kp'^ (the best rescaling factor is 1/40000). tresp decreases as K^'^ increases (being given a 
single Kf^ per K^'^). However, the controller is limited to an optimum tresp of about 2 ms 
as shown by the plateau reached for ^ 10~^ [3] (arrow in figJT]). 

So far, the origin of the saturation remains unclear. Nevertheless, a brief analysis of 
the response function of the controller to a step wherein the contribution of the RMS-to- 
DC converter is neglected (c/. appendix, section B) emphasizes that the dynamic behavior 
can reasonably be predicted (triangles in figJT]) up to 2 ms. The best agreement between 
the experimental results and the model is found when considering the weak overcritically 
damped regime, namely : 

1 

^^^^^ ^ c+^c^-^K,kW ^^^^ 



with 



c = X + ^^^'^v ) (33) 



The origin of the saturation might thus be attributed to the contribution of the RMS-to-DC 
converter. 

As expected, the shortest APIC response time is approximately 6 times longer than the 
optimal PLL locking time, ~ 0.35 ms. Thus the PLL should track frequency changes much 
faster than amplitude changes. Therefore, with PLL gains insuring a locking time much 
shorter than the APIC one, the two blocks can be considered as operating separately. Then, 
no amplitude changes which would be the consequence of a bad tracking of the resonance 
frequency can occur. 

It might be objected that the experiments and the analysis, despite consistent, have been 
performed without considering the tip-sample interaction. Regarding the PLL, the way the 
dynamics is affected when the tip is close to the surface has not yet been investigated. But 
regarding the amplitude controller. Couturier et al. [s] have reported a theoretical analysis 
of the controller stability upon the gains and the strength of the non-linear interaction in 
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the self-excitation scheme. The analysis stresses that the stability domain of the controller 
shrinks when the contribution of the non-linear interaction (pure Van der Waals) increases. 
Thus, a couple {K^\Kf^) initially inside the stability domain might correspond to an un- 
stable behavior of the controller close to the surface, thus introducing apparent dissipation^. 
Nevertheless, considering their parameters with a tip-surface distance ranging from infinity 
down to 0.8 A (corresponding to A/ < —250 Hz), that is very close to the surface for 
operating in nc-AFM^-, the stability domain weakly shrinks^. A similar analysis for the 
PLL-excitation scheme is still lacking and should be performed for quantitative comparison 
and discussion. But comparing their analysis to the tip-surface distances and frequency 
shifts which are being used in this work, we believe that the contribution of the non-linear 
interaction to the APIC dynamics remains weak and thus would not change drastically its 
coupling to the PLL. This point is strengthened by the results given in the following section 
dVCT]) . 

C. Apparent dissipation 

1. Contribution of the PLL gains 

Section IVB II has proved that the PLL gains were controlling the PLL locking time. 
Within the locking range, the higher K^Kd-, the faster the PLL. The test performed here is 
to compute approach curves for various values of K^Kd- Except K^Kd, the parameters are 
similar to those given in figlSl In particular, Q = 30000, K^'^ = 10~^ and K^'^ = 10"'^ s~^, 
corresponding to t^esp = 2 ms. 4 sets of K^Kd have been used, namely : 11000, 5000, 1000 
and 100 rad.s"^, corresponding to locking times of ~ 0.2 ms, 0.35 ms, 1.8 ms and > 4 ms, 
respectively. Note that the two later values are almost similar or larger, respectively, than 
tresp- The 3 kHz FIR low pass filter has been used and the frequency tracker has been 
engaged. The results are reported in figlH 

With the four sets of data, no effect on the frequency shift is observed. For K^Kd = 11000 
and 5000 rad.s"^, changes in phase, amplitude and damping are noticeably similar. The 
phase and the amphtude remain constant and subsequently, no damping occurs. On the 
opposite, for K^Kd = 1000 and 100 rad.s"^, that is for a PLL locking time of about or larger 
than the APIC one, the changes are more pronounced. With KoKd = 1000 rad.s"^ (set 3), 
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the phase strongly varies along the repulsive part of the interaction, from —93 to —79 degrees. 
Consequently, the amplitude starts dropping and the damping increases. This trend is more 
pronounced with K^K^ = 100 rad.s"^, for which the phase in the repulsive regime reaches 
—60 degrees, requiring the APIC to produce ~ 14% more excitation. According to the 
discussion put forward in section HVl such an effect is expected to be detected if it would 
occur. 

Therefore, if the PLL does not lock the incoming signal fast enough, that is for locking 
times of about or larger than 1 ms, corresponding approximately to tresp/S, an undesirable 
phase shifted signal is produced, resulting in an amplitude decrease and producing significant 
apparent dissipation. 

It's peculiar to notice that, when the cantilever is driven out of resonance as shown 
with the phase changes with K^Kd = 100 rad.s"^, no abnormal frequency shift occurs. 
As a matter of fact, close to the resonance, the phase changes of the free cantilever scale 
as, to first order : (p 7r/2 — 2Q{u — 1), where u = f/fo- Considering /o = 150 kHz 
and 6ip = ±30 degrees, then 6f = fo6ip/2Q = ±1.3 Hz, which is not visible in figIH](a). 
Similar effects have been reported by H.Holscher et al. [52]. This effect should be more 
(less) pronounced with low (high) Q values (±8 Hz or ±0.08 Hz with Q = 5000 or 500000, 
respectively) and the apparent dissipation higher (lower). 

2. Contribution of the frequency tracker 

As mentioned in paragraph HHP 5t the frequency tracker updates Ucent with the goal to 
prevent the phase due to the frequency shift, ippw, be added to the NCO output. Figure 
[9] reports two approach curves computed upon the frequency tracker is engaged or not. 
When it is disengaged, the phase continuously decreases along the attractive part of the 
interaction potential, as expected from equl20l At a tip-sample separation corresponding 
to the minimum of the interaction potential, r = 2.35 A, the phase reaches —120 degrees, 
meaning that the oscillator is then seriously driven out of resonance. Following the phase 
change, the amplitude continuously decreases and the damping strives to compensate the 
amplitude reduction, thus reaching 15% of the intrinsic damping of the free oscillator. Here 
again, such an effect should be measurable. When the tip is further approached towards 
the surface, the repulsive regime makes the frequency shift increasing and so does ippw. The 
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amplitude increases back to reach A^^ and the damping is obviously reduced. An amplitude 
growth is not expected in the repulsive region of the interaction potential, but it is the 
consequence of the bad tracking of the resonance frequency. On the other hand, when 
the tracker is engaged, as already mentioned, the phase remains constant and no apparent 
dissipation occurs. 

D. Scan lines 

In addition to the analysis of the time constant of the various blocks, it is important to 
focus at variables changes when the tip is scanned along a surface. Two types of surfaces 
have been investigated : 1- a sinusoidally corrugated surface with a spatial wavelength of 
6.6 A and a corrugation of ±0.1 nm, consistently with the lattice constant of KBr, a sample 
regularly used in the group, and 2- a surface with two opposite steps with a step height of 
3.3 A. The steps are built out of arctan functions and spread out laterally over 5 A. The 
upper terrace spreads out over 3 nm (c/. insets in fig lT^ . 

The results shown here have all been obtained by A/ regulation. The scan lines have 
been initiated from the approach curve shown in fig|3]with A/set = —60 Hz, corresponding 
to an initial tip-surface separation of about 5 A, that is in the short-range regime (c/. 
figlS^a))- The gains of the distance controller have been chosen in order to insure a critically 
damped response of the controller to a step of —1 Hz when A/set is reached, namely K^'^ = 
2 X 10~^ nm.Hz"^ and Kf'^ = 2 nm.Hz~^.s~^. For all of the following curves, the frequency 
tracker has been engaged. Three sets of parameters have been varied : the PLL gains, the 
scan speed and the APIC gains. 

1. Contribution of the PLL gains 

Paragraph IV B II has proven how and Kd gains were controlling the PLL locking 
time. Figure [TO] shows scan lines computed for three values of KqK^, namely : 100, 1000 
and 5000 rad.s"^, corresponding to locking times > 4, 1.8 and 0.35 ms, respectively. The 
unchanged parameters are : scan speed = 7 nm.s~^, K^'^ = 10^^ and K^'^ = 10^^ s^^. The 
latter gains correspond to tresp — 2 ms (arrow in fig|7]). Each signal, namely topography, 
A/, ip, Arms and relative damping, consists of 256 samples. The topography signal follows 
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accurately the surface corrugation. No contribution due to the gains is revealed. A/ is 
modulated around A/set = —60 Hz, with an amplitude ranging from ±2 to ±3 Hz upon 
KqK^. The accuracy of the distance control is then of about 97%. Note also that the 
nonlinear interaction makes the modulation asymmetric around A/set and the maxima are 
mismatched compared to the maxima of the surface. However, the mismatch does not 
depend on K^Kd- If the PLL is slow, a rather important phase lag is observed, ranging from 
—65 to —110 degrees, resulting in small amplitude changes. Here also, the asymmetry around 
—90 degrees is manifest and it's interesting, despite expected, to notice the doubling in the 
periodicity of the amplitude fluctuation. The related relative apparent damping fluctuates 
accordingly, reaching about 12%ro, which should be experimentally detectable. 

2. Contribution of the scan speed 

Five scan speeds ranging from 1 to 20 nm.s~^ have been used, accordingly to typical exper- 
imental values for such a scan size. The unchanged parameters are : K^Kd = 5000 rad.s"^, 
= 10^^ and Kf'^ = 10~^ s~^. The surface with opposite steps has been used and each 
signal consists of 256 samples. The results are given in figJTTl At high speed, the topography 
channel starts being distorted as a consequence of a bad distance regulation as shown by the 
large A/ variations. The phase varies accordingly, but within a narrower domain. Therefore 
neither relevant variations of amplitude nor of relative damping (±2% only) are revealed. 

3. Contribution of the A PIC gains 

Similar experiments have been carried out by varying the APIC gains. Thirteen sets of 
Kp'^ and K'^'^ values have been used over 2 orders of magnitude for each gain. The unchanged 
parameters are : scan speed = 5 nm.s~^ and K^K^ = 5000 rad.s^^. The surface with opposite 
steps has been used. For those curves, each signal consists of 1024 samples. The sets of gains 
have been chosen such that the response time of the controller is varied from 2 to 20 ms, 
according to figE The results are reported in figHJl Here again, the topography signal 
accurately follows the corrugation. In particular, no unwanted overshoot is observed at the 
step edge despite A/ varies significantly. A small phase variation is observed at the step 
edge. For the latter three channels, no dependence is observed upon the gains such that the 
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curves match each other. The small phase variation induces a tiny amplitude change, barely 
visible in the inset of fig JTW d). but the overall fluctuations are weak, which corresponds to 
relative damping fluctuations of ±3% (worse case, t^^sp ~ 20 ms). Following discussion of 
section [TV] this value is below the threshold limit of thermal noise, at least for experiments 
carried out in UHV and at room temperature. Upon gains, a small spatial shift is observed 
in the amplitude or in the relative damping signals up to a maximum value of 0.1 A. 

VI. DISCUSSION 

The above analysis stresses flve important results : 

• The PLL dynamics plays a major role in the occurrence of relevant apparent damping 
if the locking time is about or larger than 1 ms, that is only twice faster than the 
APIC optimum response time. By "relevant damping" , we mean, on the base of the 
discussion given in section IIVI a damping which would be detectable experimentally 
upon operating conditions (c/. table [T]). 

• The frequency tracker, the aim of which is to update the PLL center frequency to make 
it matching the actual resonance frequency, plays also a major role in the occurrence 
of apparent damping. It has to be mandatorily engaged when performing approach- 
curves, otherwise unwanted additional phase shift due to the PLL occurs and the 
cantilever is then driven off resonance. 

• The PLL optimal locking time is about 0.35 ms that is 6 times shorter than the 
shortest APIC response time of the free cantilever. Therefore the resonance condition 
is expected to be always properly maintained. Consequently, when the PLL operates 
properly, no amplitude changes due to a bad tracking of the resonance frequency are 
expected to occur. If they would, this should rather be the consequence of the APIC 
and/or the DPIC dynamics. 

• The APIC response time seems to be limited to tresp — 2 ms due to the RMS-to-DC 
converter. There is a priori no fundamental restriction to the APIC response time, 
as shown by flgE However, it must be stressed that if the APIC is made faster, the 
PLL should be made faster accordingly. 
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• A weak contribution of the APIC to apparent dissipation is observed. Although spatial 
shift and apparent dissipation can conditionally be generated, the overall strength of 
the effect remains weak and should hardly be measurable for UHV investigations at 
room temperature. 

In order to compare our results to other works, the contribution of the noise (equ l30l) 
to the dissipation has been estimated with the parameters given by Gauthier et al. 
The maximum of relative apparent damping they report is about 5% (c/. fig. 4 in the 
above reference, curve tl6), corresponding to a scanning speed of about 90 A.s~^. The 
spatial wavelength of their surface model being SA, the bandwidth of the corresponding 
measurement is -B ~ 11 Hz (c/. equ JM]) . The non-linear interaction is depicted by a 
Rydberg function with Uq = A x 10"" J, A = 0.599 A and = 3 A. With Aq = 1.5 nm, 
/o = 159154 Hz, kc = 26 N.m'^, Q = 24000 (E^^ = 4.8 x lO'^ eV/cycle), we get 5Ed/Ed, ~ 
2.3% at a distance^ s = 4 A. Therefore, due to the use of a low amplitude^, the strength 
of the effect would not be balanced by the noise and could be detectable experimentally. 
On the opposite, considering smaller scan speeds, the maximum of apparent dissipation 
decreases down to 2%. Such effects become then unlikely to be observed experimentally. 

Finally, let's note that the contribution of the third controller, the DPIC, has not been 
assessed in this work. Nevertheless, as mentioned by H. Hug and A. Baratoff [31], in or- 
der to minimize feedback errors and resulting image distortions, the time constant of the 
distance controller must be shorter than the speed in the fast scanning direction / lateral 
extent of the smallest feature to be resolved and necessarily (much) larger than other time 
constants (RMS-to-DC, APIC, PLL). If those conditions are satisfied, then its contribution 
to the overall stability of the setup should be weak. This is what is readily seen in figsJTOl 
[T2] and, to a certain extend in figJTTl with reasonable speeds. This is also confirmed by 
Couturier et al. who have recently found out that the distance controller had no effect on 
the stability diagram in the self-excitation scheme and that only the amplitude controller 
played a role^*^. 

To summarize, this analysis has emphasized that a maximum of about 15% of apparent 
dissipation, mainly due to the PLL and not to the APIC, could be generated (c/. figJTOl). 
The contribution of the APIC, within the range of gains used, is systematically smaller. We 
finally infer that, with our setup (UHV and room temperature) and under typical experimen- 
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tal conditions (scanning speed, PLL and APIC gains...) corresponding to Ea^ — 2 eV/cycle, 
apparent dissipation in the range of a few eV per cycle, as frequently reported in the litera- 
ture, is unlikely to occur. A striking example is provided by R. Hoffmann et al. [l^ who pu t 
in evidence a dissipation of 3 eV/cycle on a NiO(OOl) sample, despite operating at 4 K 57|. 
This strongly suggests the interpretation of the experimental damping images in terms of 
"physical" dissipation. 



VII. CONCLUSION 



The realization and successful testing of a modular nc-AFM simulator implemented with 
LabView"""^ is reported. The design is based on a real electronics which includes a digital 
PLL, the output of which is used to detect the frequency shift but also to generate the time 
dependent phase of the excitation signal. Good agreement is obtained between the locking 
behavior of the real PLL and the PLL from the simulator. The optimum locking time of the 
PLL is found to be about 0.35 ms. The behavior of the amplitude controller is also found 
to correctly describe the real setup with an optimum response of 2 ms. The analysis of the 
time constants of the former two components provides evidence that the electronics tracks 
properly the cantilever dynamics if the PLL runs more than twice faster than the amplitude 
controller. When the system is operated with properly chosen parameters, frequency shift 
vs. distance curves successfully compare to an analytic expression which ignores the finite 
response of the electronics. No phase deviation resulting in apparent dissipation occur if 
the center frequency of the PLL tracks the resonance frequency shifted by the tip-sample 
interaction. An estimate of the minimum dissipation expected to be detected experimen- 
tally gives some insights on the relevance of apparent dissipation which can conditionally 
be generated numerically. This provides a framework to discuss the overall contribution of 
apparent dissipation during experiments. Computations of scan lines show that when the 
system is operated with experimentally relevant parameters, the contribution of the propor- 
tional and integral gains of the amplitude controller and the scan speed (up to 20 nm.s"^) 
do not lead to significant apparent dissipation. To give orders of magnitude, the worse situ- 
ation (frequency tracker disengaged) leads to a maximum of 15% more dissipation than the 
intrinsic dissipation of the free cantilever. This is below the values which are experimen- 
tally reported. This strongly suggests that nc-AFM damping images mainly reflect physical 
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channels of dissipation and not electronics artifacts. 

Besides Ref.— , two publications dealing with the implementation and/or performance of 
"virtual force microscopes" have appeared recently. These simulation codes are analogous 
to ours but differ in detail. Kokavecz et al.^*^ proposed and tested a numerical scheme de- 
signed to produce response times of the whole simulated setup, as well as of separate blocks 
(amplitude, phase and distance controllers) as short as 0.1 ms. Trevethan et al.— used the 
scheme described in Ref.'^ to compute fingerprint-like responses in the frequency shift, the 
minimum tip-sample distance and the damping signal caused by an atomic-scale configu- 
ration change at the surface of the sample. This change was first predicted from atomistic 
simulations and then induced upon approach under distance control down to a judiciously 
chosen frequency setpoint. A manual describing the combined atomistic simulation^fi and 
virtual force microscope codes is now available online^. 
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Appendices 



A. The analytical method 



The analytical method gives analytical, tractable expressions of the frequency shift versus 
distance upon the force expression. Couple of years ago, F.J. Giessibl and M.Guggisberg 
have proposed approached expressions for the frequency induced by a Morse potential^i^. 
Here a complete expression is provided. The following calculation is based on a variational 
method extensively detailed in former articles^S. We start from equation [2] and postulate a 
solution of the differential equation where the amplitude and the phase are not constant, 
but are allowed to vary slowly within time, that is over durations much longer than the 
oscillating period, namely : 



z{t) = A{t) cos [ut + ip{t)] 



(34) 



Equation [2] is then equivalent to 



with : 



a cos {'^(t)) — j3 sin {(p{t)) = 7 
f3 cos {ip{t)) +asm{(p{t)) = 



(35) 



(36) 



a = A{t) + ^i(t) + {uj^, -[UJ + m?} Mt) + 
(3 = Ait)m + [2i(t) + ^A(t)] [cu + ^(t)] + 

Fc and Fd are the first Fourier components of the conservative and dissipative forces expe 
rienced by the tip over one oscillation period T, respectively : 



1 
r 

1 



T 



F, = -dA^j^ VUt)dt\ (37) 
djrVd{t)dt} (38) 



Setting A{t) = A(t) = (pit) = (p{t) = yields to the steady equations of the oscillator in 
amplitude and phase out of which the relevant variables, namely A/ = / — /o and driving 
amplitude, related to the damping, can be extracted : 
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F-kA(^-l^^^] (39) 

2/0 \^A(l + f)(l + l/tan(v.)) 

Equation [39] illustrates that, if the phase is properly maintained at —tt/2 rad, the frequency 
shift is essentially due to the conservative interaction, whatever the dissipative contribution. 
Besides, equj39]gives the well known expression of the frequency shift due to the conservative 
contribution^ : 

^ = ^ (41) 

The expression of Fc related to the Morse potential given in equE] can be integrated analyt- 
ically : 

= / " V^-^m = (2To,vA - To,2/a) , (42) 

Jo ^ 

where : 



T„^^ = e-'^(^-'^^) X Bessell(a, (3A) (43) 

Bessell(a, (3A) is the modified Bessel function of first kind for the parameters a and (3 A. As 
shown by equj37l the method not only requires to estimate the action related to Sf^"^^"^, but 
also its derivative with respect to the amplitude A : 

dAST- = -^{T,,,/,-T,,y,) (44) 
Regarding the Van der Waals contribution, it was demonstrated^ : 

Q r-vdw _ nHRA 

~ 3.iD^-A^f' ^''^ 
The frequency shift due to both contributions can finally be deduced : 



fo KA 



HRA 2U( 

6 (D2 - A2) 



372 + ^ ('^l.l/A - '^1,2a) 



(46) 
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B. APIC response time 

The problem consists in analyzing the dynamics of the oscillator the amplitude of which 
is controlled by a proportional-integral controller. For that purpose, all the components are 
assumed to be linear. No interaction occurs between the oscillator and the surface. The 
contributions of the RMS-to-DC and of the band pass filter are neglected. The system under 
investigation is given in figure [131 



1. Transfer function of the closed loop 

We start from the set of equations [35] wherein the phase is assumed to be constant and 
fixed to ip{t) — > = — 7r/2. This assumption is argued by the PLL behavior which maintains 
the phase almost constant, even when the interaction occurs (c/. fig. [3](b)). The assumption 
implies : u = uq, ip{t) = ip{t) = Q and obviously, we only focus at changes occurring in the 
resonance amplitude A{t) = Ao(t). The set of equations [33 is then equivalent to : 

= 7 

(47) 

a = Q 

Keeping = 7 yields to : 

io(t) = ^KA_-^j (48) 



2 1^ Q 

The transfer function of the block standing for the oscillator can thus be written : 
where : 



b = ^ 
" 2 

" 2Q 



(50) 



The transfer function of the APIC being Gapic{s) = K^^ + / the transfer function of 
the closed loop Gcl{s) = KiAq{s) / Aq*^-^ (s) can now be calculated : 

^ , . (51) 

s^ + ia + bK^'^)s + bKf^ 
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with 



= K,K,K- ^^^^ 

The proportional and integral gains are scaled by the transfer functions of the piezoelectric 
actuator and of the photodiodes, and Ki, respectively. 

2. Analogy 

Equation EH has two poles : 



si^2 = -c±^c^-bKf', (53) 
where c is the parameter given in equation [331 which can also be written : 

a + bk^ , , 

c = ^ (54) 

Equation [5l] is thus almost analog to a standard 2"^^ order system, the canonical form of 
which can be written : 

^(^) = 2^nf'' ^ 2 ^ (55) 

where ujn and C, are the characteristic frequency and damping factor of the system, respec- 
tively. Thus : 



C = ^= (56) 

'bKr 



Now, it's well known that the position of Q with respect to 1 defines the overall behavior of 
the system : 



Undercritically damped regime =^ C, < 1 <^ c < \/ bKf 



Critically damped regime ^( = l^c= \ bk^ (57) 



Overcritically damped regime ^ C > 1 -v^ c > ybKf 
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3. Analysis to a step response 



To assess how fast the controller reacts, we investigate the response of the controller 
to a step in amplitude upon it is in the over-, under- or critically damped regime. Let's 
assume a step of amplitude As, the corresponding transfer function is Gs{s) — Ag/s and 
the transfer function of the closed loop system to which the step is applied is therefore 
Gds{s) = Gs{s) X Gci{s), that is : 



Gcis{s) 



AsbK^ 



(s + c) 



c2 + bKf^ 



+ 



A.bKr' 



X 



s + c) - c2 + s 



The time- dependent solutions Qdsii) = C ^{Gcis{s)} are : 



• Overcritically damped regime : c> y hK'i 



(58) 



with ^± = (c ± C - hK^) I2i and C = a/c2 - hKf. 



Critically damped regime : c— \j bKf^ 

gdsit) ^ As [l - e-^' + 



bK^ - c 



te 



-ct 



(59) 



(60) 



Undercritically damped regime : c < ^bKf 



gcis{t) =As<l-e ^ X 



cos(e^t)+ " sin(rt) 



with ^' = \ bKf^ - c2. 



(61) 



4- Summary 

The transition from the over to the under critically damped regime is thus controlled by 
the parameter c, that is the proportional gain, and occurs when the condition c = \JbRf^ 
is fulfilled. In the critically damped regime, the relationship between the two gains is : 
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The time constant of the system can be extracted upon the regime and the time scale that 
are considered. For short time scales, the response time of the controller is typically 1/c in 
the critically damped regime and tresp = V + — bK^'^) in the overcritically damped 
regime (c/. equ J59|) . which is given in equation [321 Figure [7] shows that the response times of 
the real system and of the simulated machine reasonably match tresp while the RMS-to-DC 
converter has a negligible influence on the system dynamics. 
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Figures 



PLL 
board 




FIG. 1: Scheme of the simulator operating in nc-AFM, based on the design of the electronics of 
the real apparatus. 
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FIG. 2: Phase shift between a 150 kHz sinusoidal waveform sent to the real PLL and the PLL 
output waveform (/cent = 150 kHz) when tuning the input frequency from —150 to +150 Hz upon 
the frequency tracker is engaged or not. The 3kHz FIR low pass filter has been used. When the 
tracker is disengaged, the phase lag can reach ±80 degrees (continuous black line, amplitude of the 
input waveform : = 110 mV peak-to-peak). When it is engaged, the phase lag drops to almost 
zero (dotted black line and inset, similar input waveform). When the input frequency accurately 
matches the center frequency, the shift is zero (modulo an error corresponding to a few fractions 
of degrees as shown in the inset). 
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FIG. 3: Numerical approach curves. The parameters are Af^^ = 7 nm, /q = 150 kHz, k^. = 
30 N.m-i, Q = 30000, therefore Tq = 31.4 s'^ and E^o = 0.96 eV/cycle, Ki = 0.1 V.nm"!, 



K3 = 1/Ki nm.V-\ Kd = 1 V, iT, 



^1 0-1 



10" 



Kf^ = 10 s ^, approach 



m ■ 



speed 2 nm.s ^. The parameters of the interaction potential have been taken from ref. 
H = 1.865 X 10"^^ J, R = 5 nm, Uq = 3.641 x 10"^^ J, A = 1.2 A, and = 2.357 A. Except in 

(a) , the signals are monitored at 10 kHz. (a)- Comparison between the A/ computed numerically 
(open circles) and the analytic expression of A/ (thick grey line) due to Van der Waals and Morse 
interactions (equ]46p. The two curves match accurately along the attractive and repulsive parts 
of the interaction potential. For clarity reasons, 10 times less samples are displayed compared to 
plots shown in (b), (c) and (d). The dotted line depicts the analytic A/ due to a pure Van der 
Waals potential, thus showing where the short- and long-range interaction regimes are discernable. 

(b) - Phase lag, (p. The frequency tracker being engaged, the phase remains constant and equal to 
—90 degrees within deviations limited to 0.3%, thus maintaining the cantilever driven on resonance. 

(c) - Amplitude Aj-^sit)/Ki. Since no phase variation occur, the amplitude remains constant as 
well throughout the approach, (d)- Relative damping. 
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FIG. 4: (Color online) Step response of the real (a) and of the simulated (b) PLLs to a center 
frequency step of +10 Hz at /cent = 150 kHz, resulting in a A/ of —10 Hz. For this experiment, no 
interaction between the tip and the surface occurs. The various curves represent the experiments 
carried out for various values of the related gains of the PLLs denoted by the symbols. The PLL 
output is recorded at 10 kHz. The curves are fitted with a decaying exponential (thick continuous 
lines) out of which the PLL "locking time" is extracted. They are displayed over similar relative 
ranges. 
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FIG. 5: Locking time of the simulated (filled squares) and real (empty circles) PLL vs. K^Kd and 
loop gain X 91000. The locking times are obtained from the related step response curves (figsjJKb) 
and (a), respectively). The error bars depict the uncertainty on the fitted value of the locking time 
(±10%). The arrow indicates the value of the loop gain used experimentally which corresponds to 
an optimum behavior of the PLL and a related locking time of about 0.35 ms. The curve is given 
as guide eyes. 
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Time (s) Time (s) 

FIG. 6: Step response of the real (a) and of the simulated (b) APIC to a ^q''* step. To perform 
the experiments, the related PLLs are engaged (frequency trackers as well). The 3 kHz FIR low 
pass filter has been used. No interaction between the tip and the surface occurs. The cantilever 
properties are /o = 157514.6 Hz and Q = 36000. To perform the calculation, since the cantilever 
stiffness was not accurately known, we have arbitrarily chosen kc = 30 N.m^^, in reasonable 
agreement with manufacturer's datasheet. The other numerical parameters are similar to those 
given in figEl The curves depict the experiments carried out for various values of K^'^ and Kf^ 
gains. Three behaviors are observed : overcritically damped responses without overshoot, critically 
damped responses with a slight overshoot and undercritically damped responses with an oscillating 
behavior. The insets show them for couple of curves which have been arbitrarily shifted along the 
time axis, but note that the relative ranges are similar. 
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FIG. 7: Response time of the APIC vs. Kp'^ of the simulated setup and the rescaled Kp gain of the 
real controller. The best agreement between the curves is obtained with i^p/40000. The two curves 
match with a reasonable agreement and exhibit two domains : first the response time decreases 
when increasing and then a saturation is reached corresponding to tresp — 2 ms. The dotted 
line is given as guide eyes. Such an analysis can be performed assuming that the step response is 
governed by a single time constant, thus restricting the analysis to curves which exhibit an almost 
critically damped behavior (c/. text). The triangles depict the trace of the function tresp (equ. 
32]) with /o = 157514.6 Hz, Q = 36000, Ki = 0.1 V.nm^^ and K3 = 1/Ki nm.V"^ The arrow 
denotes the value of the gain used to perform the scan lines (c/. section |V]). Beyond Kp^ = 10^^, a 
noticeable discrepancy between the theoretical model and the experiments is observed which might 
be due to the contribution of the RMS-to-DC converter. 
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FIG. 8: Approach vs. distance curves for various sets (1 to 4) of K^K^ of the Ph gains, namely : 
11000 (continuous black line), 5000 (continuous grey line), 1000 (dotted black line) and 100 rad.s"^ 
(dotted grey line), corresponding to locking times of 0.2, 0.35, 1.8 and > 4 ms, respectively, (a)- 
A/, the curves are all matching each other, (b)- Phase lag tp, (c)- amplitude Arms{t) / Ki and (d)- 
relative damping. Except KoKd, the parameters are similar to those given in figlU in particular 
Kp^ = and Kf^ = 10^^ s^^ corresponding to trosp — 2 ms. When the PLL locking time is 
larger than Uesp (set 4, KqK^i = 100), the resonance is not properly locked, which induces a phase 
shift. Consequently, the amplitude decreases and the damping increases. 
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FIG. 9: Approach vs. distance curves upon the frequency tracker of the PLL is engaged or not (grey 
or black hues, respectively). (a)- Frequency shift A/, (b)- phase lag ip, (c)- amplitude Aj-^s(t)/Ki 
and (d)- relative damping. The parameters are similar to those given in the caption of figl3l When 
not engaged, the phase continuously drifts during the approach due to the increase of the attractive 
interaction meaning that the oscillator is not driven on resonance. Subsequently, the amplitude 
drops and the APIC strives to keep it constant by increasing the excitation. 15% more excitation 
is thus produced, that is above thermal noise (c/. section HV]) . 
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FIG. 10: Calculated cross-section of a sinusoidally corrugated surface for various PLL gains, namely 
KoKfi = 100 rad.s"^ (continuous black line), 1000 rad.s"^ (continuous grey line) and 5000 rad.s"^ 
(dotted black line). The scan lines have been initiated from the approach curve shown on figl3]by 
A/ regulation using A/set = —60 Hz. The lateral scan speed is 7 nm.s^^ and the section consists 
of 256 samples. Kp^ = 2 x 10^^ nm.Hz^^ and Kf"^ = 2 nm.Hz^^.s^^, corresponding to a critically 
damped response of the controller to a frequency step of —1 Hz. (a)- Topography, (b)- Frequency 
shift, (c)- Phase lag if. (d)- Amplitude Arms/-f^i and (e)- relative damping. No noticeable effect is 
revealed on the topography. The apparent dissipation remains below the thermal noise (c/. section 
HVll . except if the PLL is slow {KoK^ = 100 rad.s-\ 12%). 
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FIG. 11: Calculated cross-section of a surface with two opposite steps for various scanning speeds, 
namely 1 (■), 2 (•), 5 (a), 10 (+) and 20 nm.s^"'^ (o). For clarity reasons, the right hand side 
step region has been magnified. The insets show the whole section. The scan lines have been 
initiated upon the same conditions than in fig [10] with K^K^ = 5000 rad.s^^, K^^ = 10~^ and 
Kf" = 10""^ s-^ The high speeds require to reduce the number of samples per line to 256. The 
gains of the distance controller are similar to those given in figllOl (a)- Topography, (b)- Frequency 
shift, (c)- Phase lag ip. (d)- Amplitude Arms/^'i and (e)- relative damping. At high scan speeds, 
the topography is slightly distorted, but the overall damping remains weak and below the thermal 
noise (c/. section HV|) . 
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FIG. 12: Calculated cross-section of a surface with two opposite steps for various sets of the 
APIC gains, namely {K^" = lQ-'^]Kf^ = lO^^ g-^) = o, (10-3;10-3) = o, (IQ-^jlQ-^) = x, 
(10-3; 10-2) = +, (10-2; 10-4) =A, (10-^; IQ-^) = ., (IQ-^ 10-^) = ■, (IQ-^ 10-^) = 
(10-3; 10-^) = A, (10-2; IQ-S) = ▼ and (IQ-^; IQ-^) = 0. The scan lines have been initiated 
upon the same conditions than in fig(Tn] with KqK^ = 5000 rad-s^^. The lateral scan speed is 
5 nm.s-^ and the section consists of 1024 samples. The gains of the distance controller are similar 
to those given in figllOl (a)- Topography, (b)- Frequency shift, (c)- Phase lag 99. (d)- Amplitude 
^rms/^i and (e)- relative damping. The APIC gains have a negligible effect on the topography. A 
weak apparent damping (3%ro, that is not experimentally relevant) is revealed at the step if the 
APIC is slow, corresponding to tresp = 20 ms. 
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FIG. 13: Simplified scheme of the closed loop for the characterization of the response time of the 
APIC. 
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